Simulation of Heavily Irradiated Silicon Pixel Detectors 
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We show that doubly peaked electric fields are necessary to describe grazing-angle charge collection measure- 
ments of irradiated silicon pixel sensors. A model of irradiated silicon based upon two defect levels with opposite 
charge states and the trapping of charge carriers can be tuned to produce a good description of the measured 
charge collection profiles in the fluence range from 0.5xl0 14 n eq /cm 2 to 5.9xl0 14 n eq /cm 2 . The model cor- 
rectly predicts the variation in the profiles as the temperature is changed from — 10°C to — 25°C. The measured 
charge collection profiles are inconsistent with the linearly- varying electric fields predicted by the usual descrip- 
tion based upon a uniform effective doping density. This observation calls into question the practice of using 
effective doping densities to characterize irradiated silicon. The model is now being used to calibrate pixel hit 
reconstruction algorithms for CMS. 



1. Introduction 

A silicon pixel detector [lj is currently being de- 
veloped for the CMS experiment at the CERN Large 
Hadron Collider (LHC). The detector will be a key 
component in the reconstruction of primary and sec- 
ondary vertices in the particularly harsh LHC environ- 
ment characterized by large track multiplicities and 
high radiation backgrounds. The innermost layer, lo- 
cated at only 4 cm from the beam line, is expected to 
be exposed to a 1 MeV neutron equivalent fluence of 
3 x 10 14 n eq /cm 2 per year at full luminosity. 

The response of the silicon sensors during the de- 
tector operation is of great concern. It is well under- 
stood that the intra-diode electric fields in these detec- 
tors vary linearly in depth reaching a maximum value 
at the p-n junction. The linear behavior is a conse- 
quence of a uniform space charge density, N c g , caused 
by thermally ionized impurities in the bulk material. 
It is well known that the detector characteristics are 
affected by radiation exposure, but it is generally as- 
sumed that the same picture is valid after irradiation. 
In fact, it is common to characterize the effects of irra- 
diation in terms of a varying effective uniform charge 
density. In Q we have proved that this picture does 
not provide a good description of irradiated silicon 
pixel sensors. In addition, it was shown that it is 
possible to adequately describe the charge collection 
characteristics of a heavily irradiated silicon detector 
in terms of a tuned double junction model which pro- 



duces a doubly peaked electric field profile across the 
sensor. The modeling is supported by the evidence 
of doubly peaked electric fields obtained directly from 
beam test measurements and presented in 0. The 
dependence of the modeled trap concentrations upon 
fluence was presented in and the temperature de- 
pendence of the model was presented in [5j . We sum- 
marize these results in this document. 

This paper is organized as follows: Section de- 
scribes the experimental details, Section [3] describes 
the carrier transport simulation used to interpret the 
data. The tuning of the double junction model and its 
resulting predictions are discussed in Section 0] The 
temperature dependence of the data and model are 
summarized in Section |SJ The conclusions are given 
in Section 



2. Experimental Details 

The measurements were performed in the H2 beam 
line of the CERN SPS in 2003/04 using 150-225 GeV 
pions. The beam test apparatus is described in |(J. 
A silicon beam telescope [3 consisted of four modules 
each containing two 300 /im thick single-sided silicon 
detectors with a strip pitch of 25 fim and readout pitch 
of 50 /mi. The two detectors in each module were ori- 
ented to measure horizontal and vertical impact co- 
ordinates. A pixel hybrid detector was mounted be- 
tween the second and third telescope modules on a 
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cooled rotating stage. A trigger signal was generated 
by a silicon PIN diode. The analog signals from all 
detectors were digitized in a VME-based readout sys- 
tem by two CAEN (V550) ADCs and one custom-built 
flash ADC. The entire assembly was located in an 
open-geometry 3T Helmholtz magnet that produced 
a magnetic field either parallel or orthogonal to the 
beam. The temperature of the tested sensors was 
controlled with a Peltier cooler that was capable of 
operating down to -30° C. The telescope information 
was used to reconstruct the trajectories of individual 
beam particles and to achieve a precise determination 
of the particle hit position in the pixel detector. The 
resulting intrinsic resolution of the beam telescope was 
about 1 /im. 

The prototype pixel sensors are so-called "n-in-n" 
devices: they are designed to collect charge from 
n + structures implanted into n-bulk silicon using p- 
spray isolation. All test devices were 22x32 arrays of 
125x125 /an 2 pixels that were fabricated by CiS. The 
substrate, produced by Wacker, was 285 /jm thick, 
n-doped, diffusively-oxygenated float zone silicon of 
orientation (111), resistivity 3.7 kfl-cm. and oxygen 
concentration in the order of 10 17 cm~ 3 . Individual 
sensors were diced from fully processed wafers after 
the deposition of under-bump metalization and in- 
dium bumps. A number of sensors were irradiated 
at the CERN PS with 24 GeV protons. The irradia- 
tion was performed without cooling or bias. The de- 
livered proton fluences scaled to 1 MeV neutrons by 
the hardness factor 0.62 0| were 0.5 xlO 14 n oq /cm 2 , 
2xl0 14 n oq /cm 2 and 5.9 xlO 14 n eq /cm 2 . All samples 
were annealed for three days at 30° C. In order to avoid 
reverse annealing, the sensors were stored at -20° C af- 
ter irradiation and kept at room temperature only for 
transport and bump bonding. All sensors were bump 
bonded to PSI30/AC30 readout chips which allow 
analog readout of all 704 pixel cells without zero sup- 
pression. The PSI30 settings were adjusted to provide 
a linear response to input signals ranging from zero to 
more than 30,000 electrons. 



simultaneously solve Poisson's Equation, the carrier 
continuity equations, and various charge transport 
models. A final simulation step reformatted the sim- 
ulated data into test beam format so that it could be 
processed by the test beam analysis software. 

The simulation was checked by comparing simu- 
lated data with measured data from an unirradiated 
sensor. A plot of the charge measured in a single pixel 
as a function of the horizontal and vertical track im- 
pact position for normally incident tracks is shown 
in Fig. ^ The simulation is shown as the solid his- 
togram and the test beam data are shown as solid 
points. Note that the sensor simulation does not in- 
clude the "punch-through" structure on the n + im- 
plants which is used to provide a high resistance con- 
nection to ground and to provide the possibility of on- 
wafer IV measurements. There is reduced charge col- 
lection from this portion of the implant and the data 
shows reduced signal in both projections at the bias 
dot. Another check, shown in Table [IJ is the compar- 
ison of the average Lorentz angle measured at several 
bias voltages |(|. In both cases, reasonable agreement 
is observed between measured and simulated data. 
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Figure 1: Collected charge measured in a single pixel as 
a function of the horizontal (left) and vertical (right) 
track impact position for tracks that are normally 
incident on an unirradiated sensor. The simulation is 
shown as a solid histogram and the test beam data are 
shown as solid dots. 



3. Sensor simulation 

The interpretation of the test beam data relies upon 
a detailed sensor simulation that includes the model- 
ing of irradiation effects in silicon. The simulation, 
pixelav 0, El EI, incorporates the following ele- 
ments: an accurate model of charge deposition by 
primary hadronic tracks (in particular to model delta 
rays); a realistic 3-D intra-pixel electric field map; 
an established model of charge drift physics including 
mobilities, Hall Effect, and 3-D diffusion; a simula- 
tion of charge trapping and the signal induced from 
trapped charge; and a simulation of electronic noise, 
response, and threshold effects. The intra-pixel elec- 
tric field map was generated using tcad 9.0 0] to 



Table I Measured and simulated values of average 
Lorentz angle 6l versus bias voltage for an unirradiated 
sensor. 



Bias Voltaj; 


;e Measured 9l [deg 


Simulated 0l [deg] 


150V 


22.8 ± 0.7 


24.7±0.9 


300V 


14.7 ±0.5 


17.4±0.9 


450V 


11.2 ± 0.5 


12.0±0.9 



The effect of irradiation was implemented in the 
TCAD simulation by including two defect levels in the 
forbidden silicon bandgap with opposite charge states 
and trapping of charge carriers. The model, similar to 
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one proposed in |l3|. is based on the Shockley- Read- 
Hall (SRH) statistics and produces an effective space 
charge density p e g from the trapping of free carriers 
in the leakage current. The effective charge density 
is related to the occupancies and densities of traps as 
follows, 



Pefi = e [N D f D - N A f A ] + Pdopants 



(1) 



where: No and N A are the densities of donor and ac- 
ceptor trapping states, respectively; fn and f A are the 
occupied fractions of the donor and acceptor states, 
respectively, and Pdopants is the charge density due to 
ionized dopants (describes the resistivity of the mate- 
rial before irradiation) . The donor and acceptor occu- 
pancies are related to the trap parameters by standard 
SRH expressions 



So 



fA 



v h rfp + v e o-V ni e ED / kT 



(2) 



v e afn + v h a£nie~ EA/kT 



where: v e and vt are the thermal speeds of electrons 
and holes, respectively; c,P, a® are the electron and 
hole capture cross sections for the donor trap; af-, 
are the electron and hole capture cross sections for 
the acceptor trap; n, p are the densities of free elec- 
trons and holes, respectively; rij is the intrinsic den- 
sity of carriers; Ejj, E A are the activation energies 
(relative to the mid-gap energy) of the donor and ac- 
ceptor states, respectively. Note that the single donor 
and acceptor states model the effects of many physical 
donor and acceptor states making the two-trap model 
an "effective theory" . 

The physics of the model is illustrated in Fig. [3 
The space charge density and electric field are plot- 
ted as functions of depth z for a model tuned to re- 
produce the $ = 5.9 x 10 14 n oq cm~ 2 charge collec- 
tion data at 150V bias. The SRH process produces 
electron-hole pairs more or less uniformly across the 
thickness of the sensor. As the electrons drift to the 
n+ implant, the total electron current increases as z 
decreases. The hole current similarly increases with 
increasing z. Trapping of the mobile carriers produces 
a net positive space charge density near the p + back- 
plane and a net negative space charge density near the 
n + implant. Since positive space charge density cor- 
responds to n-type doping and negative space charge 
corresponds to p-type doping, there are p-n junctions 
at both sides of the detector. The electric field in the 
sensor follows from a simultaneous solution of Pois- 
son's equation and the continuity equations. The re- 
sulting z-component of the electric field varies with an 
approximately quadratic dependence upon z having a 
minimum at the zero of the space charge density and 
maxima at both implants. A more detailed descrip- 



tion of the double junction model and its implemen- 
tation can be found in [2( . 
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Figure 2: The space charge density (solid line) and 
electric field (dashed line) at T — — 10° C as functions of 
depth in a two-trap double junction model tuned to 
reproduce the $ = 5.9 x 10 14 n oq cm~ 2 charge collection 
data at 150V bias. 



4. Model tuning and results 

Charge collection across the sensor bulk was mea- 
sured using the "grazing angle technique" As * s 
shown in Fig. [31 the surface of the test sensor is ori- 
ented by a small angle (15°) with respect to the pion 
beam. Several samples of data were collected with 
zero magnetic field and at temperature of — 10°C and 
— 25°C. The charge measured by each pixel along the 
y direction samples a different depth z in the sen- 
sor. Precise entry point information from the beam 
telescope is used to produce finely binned charge col- 
lection profiles. 



y axis 
p+ sensor backplane 




n+ pixel implant 
Bump bond 



Figure 3: The grazing angle technique for determining 
charge collection profiles. The charge measured by each 
pixel along the y direction samples a different depth z in 
the sensor. 



The charge collection profiles for a sensor irradi- 
ated to a fluence of $ = 5.9 x 10 14 n cq /cm 2 and op- 
erated at a temperature of — 10°C and bias voltages 
of 150V and 300V are presented in Fig^J The mea- 
sured profiles are shown as solid dots and the simu- 
lated profiles are shown as histograms. In order to 
investigate the applicability of the traditional picture 
of type-inverted silicon after irradiation, the simulated 
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profiles were generated with electric field maps corre- 
sponding to two different effective densities of accep- 
tor impurities. The full histograms are the simulated 
profile for N cS = 4.5 x 10 12 cm" 3 . Note that the 300V 
simulation reasonably agrees with the measured pro- 
file but the 150V simulation is far too broad. The 
dashed histograms show the result of increasing N e fi 
to 24 x 10 12 cm -3 . At this effective doping density, 
the width of the simulated peak in the 150V distribu- 
tion is close to correct but it does not reproduce the 
"tail" observed in the data at large y. The 300V sim- 
ulated distribution is far too narrow and the predicted 
charge is lower than the data (note that the profiles 
are absolutely normalized). It is clear that a simu- 
lation based upon the standard picture of a constant 
density of ionized acceptor impurities cannot repro- 
duce the measured profiles. 




500 1000 

Position (|im 



Figure 4: The measured and simulated charge collection 
profiles for a sensor at T = — 10° C irradiated to a fluence 
of $ = 5.9 x 10 14 n cq /cm 2 . The profiles measured at bias 
voltages of 150V and 300V are shown as solid dots. The 
full histograms are the simulated profiles for a constant 
effective doping N e s = 4.5 x 10 12 cm -3 of acceptor 
impurities. The dashed histograms are the simulated 
profiles for a constant effective doping 



N eS = 24 x 10 12 



The same measured profiles and those from bias 
voltages of 200V and 450V are shown in Fig. [5] They 
are compared with simulations based upon the elec- 
tric field produced by the two trap model. The model 
has six free parameters (JVc, Na, cr^ , cr^ , af, a£) 
that can be adjusted. The activation energies are 
kept fixed to the values of pj: E D = E v + 0.48 eV, 
Ea — Ec — 0.525 eV where Ey and Eq are the en- 
ergies of the valence and conduction band edges. The 
electric field map produced by each tcad run is in- 
put into pixel AV. The electron and hole trapping 
rates, T e and Th, are also inputs to pixelav and are 
treated as constrained parameters. Although they 
have been measured |l5| . they are uncertain at the 
20% level due to the fluence uncertainty and possi- 
ble annealing of the sensors. They are therefore al- 
lowed to vary by as much as ±20% from their nom- 



inal values. The donor concentration of the starting 
material is set to 1.2 x 10 12 cm -3 corresponding to a 
full depletion voltage of about 70 V for an unirradi- 
ated device. Because each model iteration took ap- 
proximately two days, it was not possible to use stan- 
dard statistical fitting techniques. The parameters of 
the double junction model were systematically varied 
and the agreement between measured and simulated 
charge collection profiles was judged subjectively. The 
"best fits" shown in this paper are probably not true 
likelihood minima and the calculation of eight param- 
eter error matrices is beyond available computational 
resources. Adequate agreement was achieved by set- 
ting the ratio of the common hole and electron cross 
sections ah/o- e to 0.25 and the ratio of the acceptor 
and donor densities Na/Nd to 0.40. There is a range 
of parameters in the Ni)-o~ e space that produces rea- 
sonable agreement with the measured profiles. The 
range is shown in Fig. as the solid line in the loga- 
rithmic space. If the donor density becomes too small 
(JVd < 20 x 10 14 cm -3 ), the 150V simulation produces 
too much signal at large z. If the donor density be- 
comes too large (No > 50xl0 17 cm -3 ), the 300V sim- 
ulation produces insufficient signal at large z. Since 
the simulated leakage current varies as /i ca k oc a e Nu, 
different points on the allowed solid contour corre- 
spond to different leakage current. Contours of con- 
stant leakage current are shown as dashed curves and 
are labeled in terms of the corresponding damage pa- 
rameter a where an = 4 x 10 -17 A/cm is the expected 
leakage current It is clear that the simulation 

can accommodate the expected leakage current which 
is smaller than the measured current by a factor of 
three. 

The electron and hole traps in the model should 
also contribute to the trapping of signal carriers. The 
contributions of these states to the effective trapping 
rates of electrons and holes are given by the following 
expressions 



T e - v e [a*N A (l - f A ) + <J?N D f D ] ~ v e a^N A 



(3) 



T h = v h \a%N D {l - f D ) + a£N A f A \ * v h a%N. 



D : 



D 



where it has been assumed that the trap occupancies 
are small. Because Na/Nd is assumed to be constant, 
contours of constant electron trapping rate are paral- 
lel to contours of constant leakage current in N]j-a e 
space. The best "fit" of the simulation to the mea- 
sured profiles reduced T e to 85% of the un-annealed 
trapping rate Tq for the nominal fluence fl5| . These 
contours are compared with the allowed contour in 
Fig. Eh>. It is clear that the simulation can accommo- 
date the measured trapping rate in the same region of 
parameter space that maximizes the leakage current. 

Figure^ also suggests a solution to the puzzle that 
the trapping rates have been shown to be unaffected 
by the presence of oxygen in the detector bulk |l5j 
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Figure 5: The measured charge collection profiles at a 
temperature of — 10° C and bias voltages of 150V, 200V, 
300V, and 450V are shown as solid dots for a fluence of 
5.9 x 10 14 n eq /cm 2 . The two-trap double junction 
simulation is shown as the solid histogram in each plot. 
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Figure 6: The allowed region in the No-Ve space for the 
best fit 5.9 x 10 14 n eq /cm 2 model is shown as the solid 
line in (a) and (b). Contours of constant leakage current 
are shown as dashed curves in (a) and are labeled in 
terms of the corresponding damage parameter a where 
Qo = 4 x 10 -17 A/cm is the expected leakage current 
flrjl . Contours of constant electron trapping rate are 
shown as dashed curves in (b) and are labeled in terms of 
the standard-annealed trapping rate To for the nominal 
fluence [lfj. 



whereas it is well-established that the space charge ef- 
fects are quite sensitive to the presence of oxygen in 
the material ^?|]- It is clear from Fig 03 that small- 
cross-section trapping states can play a large role in 
the effective charge density but a small one in the 
effective trapping rates: every point on the dj44 line 



produces 100% of the effective charge density but only 
the larger cross section points contribute substantially 
to the trapping rate. If the formation of the additional 
small-cross-section states were suppressed by oxygen, 
then p G ff could be sensitive to oxygenation whereas 
r \/h would be insensitive to oxygenation. This is an- 
other consequence of the observation that the occu- 
pancies fnlA °f the trapping states are independent 
of the scale of the cross sections in the steady state 
(see eq. |3J). The trapping of signal carriers is not a 
steady-state phenomenon and is sensitive to the scale 
of the trapping cross sections. 

The simulation describes the measured charge col- 
lection profiles well both in shape and normalization. 
The "wiggle" observed at low bias voltages is a sig- 
nature of the doubly peaked electric field shown in 
Fig. |21 The relative signal minimum near y ~ 700 fixsx 
(see Fig-EJ) corresponds to the minimum of the electric 
field z-component, E z , where both electrons and holes 
travel only short distances before trapping. This small 
separation induces only a small signal on the n + side of 
the detector. At larger values of y, E z increases caus- 
ing the electrons drift back into the minimum where 
they are likely to be trapped. However, the holes drift 
into the higher field region near the p + implant and 
are more likely to be collected. The net induced signal 
on the n + side of the detector therefore increases and 
creates the local maximum seen near y — 900 fim. 

The charge collection profiles at T = — 10°C for sen- 
sors irradiated to fluences of $ = 0.5 x 10 14 n cq /cm 2 
and <f> = 2 x 10 n eq /cm 2 and operated at several bias 
voltages are presented in Fig. EJa-c) and Fig. Efd-g), 
respectively. The measured profiles, shown as solid 
dots, are compared to the simulated profiles, shown 
as histograms. Note that the "wiggle" is present at 
low bias even at <E> = 0.5 x 10 14 n oq /cm 2 which is just 
above the "type-inversion" fluence. This suggests that 
a doubly peaked field is present even at rather small 
fluences. 

The double junction model can provide a reasonable 
description of the lower fluence charge collection pro- 
files using the parameters obtained with the fitting 
procedure shown in Table [H] We observe that the 
donor trap concentration increases more rapidly with 
fluence than does the acceptor trap concentration. 
The ratio between acceptor and donor trap concen- 
trations is 0.76 at the lowest fluence and decreases to 
0.40 at 5.9x 10 14 n cq /cm 2 . In addition, the fits exclude 
a linear dependence of the trap concentrations with 
the irradiation fluence. At $ = 5.9 xlO 14 n oq /cm 2 
the cross section ratio erh/c7 e is set to 0.25 for both 
donor and acceptor traps while at lower fluences we 
find afj 'af = 0.25 and a® /a® = 1 for the acceptor 
and donor traps, respectively. The simulated leakage 
current is approximately linear in fluence, but the ra- 
tio Na/Nd is clearly not constant. This may be a 
consequence of the quadratic fluence scaling of one 
or more di-vacancy states or it may reflect the fact 
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Figure 8: The z-component of the simulated electric field 
at T — — 10° C resulting from the model best fit is shown 
as a function of z for a sensor irradiated to fluences of 
0.5 x 10 14 n cq /cm 2 (a) and 2 x 10 14 n cq /cm 2 (b). 



that the two trap model with the particular choice of 
activation energies does not accurately model the de- 
pendence of the trap occupancies on leakage current. 
The allowed Nu-a e parameter spaces for the lower flu- 
ence models are much more constrained than in the 
<£>=5.9xl0 14 n cq /cm 2 case and predict the expected 
leakage current. The electron and hole trapping rates, 
r e and Th are found to scale more or less linearly with 
flucnce. 



Table II Double trap model parameters extracted from 
the fit to the data. 



$ [10 14 n oq cm- 2 ] 


0.5 


2.0 


5.9 


N A [10 14 cm" 3 ] 


1.9 


6.8 


16 


N D [10 14 cm" 3 ] 


2.5 


10 


40 


at ID [10- 15 cm 2 ] 


6.60 


6.60 


6.60 


at [10- 15 cm 2 ] 


1.65 


1.65 


1.65 


a» [10- 15 cm 2 ] 


6.60 


6.60 


1.65 


T e [KT 2 ns~ 2 ] 


2.7 


9.6 


28. 


Y h [lO" 2 ns" 2 ] 


3.6 


13. 


38. 



and $ = 2 x 10 14 n oq /cm 2 . At $ = 5.9 x 10 14 n cq /cm 2 , 
the field profile has a minimum near the midplane of 
the detector which shifts toward the p+ implant at 
lower fluences. The field has maxima at the detector 
implants as discussed in Section |3J Figure 02a) shows 
that a doubly peaked electric field is necessary to de- 
scribe the measured charge collection profiles even at 
the lowest measured fluence which is just beyond the 
"type inversion point" . The dependence of the space 
charge density upon the z coordinate is shown in Fig- 
ures |21 and Before irradiation the sensor is charac- 
terized by a constant and positive space charge den- 
sity of 1.2 x 10 12 cm -3 across the sensor bulk. After a 
flucnce of 0.5 x 10 14 n cq /cm 2 the device shows a nega- 
tive space charge density of about —1 x 10 12 cm~ 3 
for about 70% of its thickness, a compensated re- 
gion corresponding to the E z minimum and a posi- 
tive space charge density close to the backplane. The 
space charge density and electric field near the p+ im- 
plant increase with fluence. The space charge density 
is not linear in z due to the variation of the carrier 
drift velocities with the electric fields. 
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The z-component of the simulated electric field, E z , 
is plotted as a function of z in Fig. [21 for 4? = 5.9 x 
10 14 n oq /cm 2 and in Fig.Elfor $ = 0.5 x 10 14 n oq /cm 2 



Figure 9: The simulated space charge density at 

T = — 10° C as a function of the z coordinate for fluences 

of 0.5 x 10 14 n cq /cm 2 and 2 x 10 14 n cq /cm 2 . 
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5. Temperature dependence 

The temperature dependence of the charge collec- 
tion profiles was studied by accumulating data at 
T = — 25°C. The pixelav simulation includes tem- 
perature dependent mobilities, diffusion, and trapping 
rates. The tcad calculation of the electric field map 
is also based upon temperature dependent quantities 
including the bandgap energy and SRH lifetimes. The 
T = — 25°C charge collection profiles for the <f> = 
2.0 x 10 14 n eq / cm 2 and $ = 5.9xl0 14 n cq / cm 2 sensors 
are compared with the simulation in Fig. 1101 It is clear 
that the simulation correctly tracks the temperature- 
dependent variations in the measured profiles. 
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Figure 10: Measured (full dots) and simulated 
(histogram) charge collection profiles at T = — 25°C and 
several bias voltages for sensors irradiated to fluences of 
2.0 x 10 14 n cq /cm 2 (a-d) and of 5.9 x 10 14 n oq /cm 2 (e-g). 



The effect of temperature on the z-component of the 
simulated electric field at $ = 5.9 x 10 14 n oq /cm 2 is 
shown in Fig.^2f° r bias voltages of 150V and 300V. It 
is clear that decreasing the temperature also decreases 
the fields on the p+ side of the sensor and increases 
them on the n+ side. 



6. Conclusions 

In this paper we have shown that doubly peaked 
electric fields are necessary to describe grazing-angle 
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Figure 11: The simulated ^-component of the electric 
field as a function of the z coordinate at the 
5.9 x 10 14 n cq /cm 2 fluence for temperatures T = — 10° C 
and T = — 25° C. The field profiles are shown for bias 
voltages of 150V and 300V. 



charge collection measurements of irradiated silicon 
pixel sensors. A model of irradiated silicon based upon 
two defect levels with opposite charge states and the 
trapping of charge carriers can be tuned to produce 
a good description of the measured charge collection 
profiles in the fluence range from 0.5 x 10 14 n eq /cm 2 to 
5.9xl0 14 n oq /cm 2 . The model correctly predicts the 
variation in the profiles as the temperature is changed 
from -10°C to -25°C. 

The doubly peaked electric field profiles have max- 
ima near the implants and minima near the detector 
midplane. This corresponds to negative space charge 
density near the n + implant and and positive space 
charge density near the p + backplane. We find that 
it is necessary to decrease the ratio of acceptor con- 
centration to donor concentration as the fluence in- 
creases. This causes the electric field profile to become 
more symmetric as the fluence increases. The effect 
of decreasing the temperature has the opposite effect 
of suppressing the fields on the p+ side of the sensor 
and increasing them on the n+ side. 

The measured charge collection profiles of irradi- 
ated sensors are inconsistent with the linearly-varying 
electric fields predicted by the usual description based 
upon a uniform effective doping density. This suggests 
that the correctness and the physical significance of ef- 
fective doping densities determined from capacitance- 
voltage measurements are quite unclear. In addition, 
we remark that the notion of partly depleted silicon 
sensors after irradiation is inconsistent with the mea- 
sured charge collection profiles and with the observed 
doubly peaked electric fields. 

The charge-sharing behavior and resolution func- 
tions of many detectors are sensitive to the details of 
the internal electric field. A known response func- 
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tion is a key element of an optimal reconstruction 
procedure. The effective model described in this pa- 
per is being used to calculate detailed response func- 
tions that are being incorporated into a new hit recon- 
struction algorithm for the CMS pixel tracking system 
[Tsl IT^ | . This will permit the "calibration" of the re- 
construction algorithm to be tracked as the pixels are 
irradiated during LHC operation. 
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